Exploring high-density baryonic matter: Maximum freeze-out density 
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Over the past decade a striking regularity has been established in heavy-ion collisions: From 
the lowest beam energies to the highest, the yields of various hadronic species are consistent 
with the assumption of chemical equilibrium [1, 2, 3, 4, 5]. Analyses of the experimentally 
obtained hadronic yield ratios at a variety of collision energies have shown that the data can 
be well reproduced within the conceptually simple statistical model that describes an ideal 
hadron resonance gas in statistical equilibrium. Furthermore, the extracted freeze-out values 
of the temperature T and the baryon chemical potential fiB exhibit a smooth and monotonic 
dependence on the collision energy and can be simply parametrized. 

The beam energy thus plays a determining role for the thermodynamic properties of the final 
state in relativistic heavy-ion collisions. However, there is no simple relationship between the 
collision energy and the freeze-out value of the (net) baryon density: At low collision energies 
the freeze-out density increases with the energy, whereas it decreases when the collision energy is 
high due to the onset of nuclear transparency. Thus there must exist a certain range of collision 
energies within which the freeze-out values of the net baryon density displays a maximum. 

The optimal collision energy leading to this highest freeze-out density was discussed in [6] 
on the basis of the up-to-date results on the properties of the final state. It was pointed out 
there that since neither fis nor T is subject to a conservation law they may be less suitable in 
a dynamical context. Furthermore, when a first-order phase transition is present, they become 
multivalued functions of the basic mechanical variables pB (net baryon density) and £ (energy 
density) inside the mixed-phase region. It is therefore of interest to reexpress the thermodynamic 
variables in terms of those mechanical densities. Accordingly, we considered in [6] how the freeze- 
out line appears when represented in terms of the basic baryon and energy densities, rather than 
chemical potential and temperature. 

In [6], we presented the freeze-out line in terms of the net baryon density pB and the energy 
density e. We show below (in Fig. 1) the corresponding {pB,£*) representation, where the 
"excitation energy density" e* = e — itinPb is the energy density above the minimum value 
rriMPB dictated by the specified net baryon density. Thus e* has both compressional and thermal 
contributions. 
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Figure 1: The hadronic freeze-out line in the 
Pb — e* phase plane as obtained from the val- 
ues of pb and T that have been extracted 
from the experimental data [2]. The calcula- 
tion employs values of pq and ps that ensure 
(5) = and (Q) = OA{B) for each value of 
Pb- Also indicated are the beam energies (in 
GeV/N) for which the particular freeze-out 
conditions are expected at either a collider 
(red) or a fixed-target facility (blue). 
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The corresponding {pB,T) diagram, where T is the freeze-out temperature is perhaps more 
easily grasped and we show it in Fig. 2 below: 

Figure 2: The hadronic freeze-out line in the 
Pb — T phase plane as obtained from the val- 
ues of and T that have been extracted 
from the experimental data [2]. The calcula- 
tion employs values of pq and that ensure 
(S) = and (Q) = OA{B) for each value of 
Pb- Also indicated are the beam energies (in 
GeV/N) for which the particular freeze-out 
conditions are expected at either a collider 
(red) or a fixed-target facility (blue). 



These novel representations of the freeze-out line bring out very clearly that there is a 
maximum value of the net baryon density: At the highest collision energies, freeze-out occurs 
for a negligible value of pB and at an energy density of nearly one half GeV/fm^; then, in the 
range of p,B = 400 — 500 MeV (and a temperature of T = 140 — 130 MeV), the freeze-out line 
exhibits a backbend and approaches the origin. Thus, the net baryon density at freeze-out has 
a maximum value which amounts to about three quarters of the familiar nuclear saturation 
density of pQ 0.16 fm~'^. 

The fact that the freeze-out value of the net baryon density exhibits a maximum as the 
collision energy is being scanned suggests that the corresponding collision energy (range) is 
optimal for the exploration of compressed baryonic matter. This suggested optimal beam kinetic 
energy is 15 — 30GeV per nucleon for a fixed-target configuration (such as FAIR at GSI), 
corresponding to ^snn = 5.6 — 7.8 GeV for a collider (such as RHIC at BNL). 

The results presented here should provide valuable guidance for establishing the desired 
capabilities of the contemplated NICA at JINR. In particular, our results suggest that freeze- 
out densities all the way up to the maximum value could be explored at a collider facility 
delivering beam kinetic energies of up to ~2.4 Ge\ /N . 
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